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We investigate the spectral statistics of the asymmetric rotor model (triaxial rigid 
rotator). The asymmetric top is classically integrable and, according to the Berry- 
Tabor theory, its spectral statistics should be Poissonian. Surprisingly, our nu- 
merical results show that the nearest neighbor spacing distribution P{s) and the 
spectral rigidity h.^(L) do not follow Poisson statistics. 



1 Introduction 

As is well known, in the semiclassical limitEI'i there is a clear connection 
between the behavior of classical systems (regular or chaotic) and the corre- 
sponding quantal ones. In particular for quantal systems, corresponding to 
classical regular systems, the spectral statistics (P(s) and A3(L)) follow the 
Poisson ensemble, while for systems corresponding to chaotic ones the Wigner 
ensemble is followed (see, for example, Ref. 3 and references therein). 

Nevertheless, some exceptions are known. The most famous case is per- 
haps the harmonic oscillator one, discussed in great detail in Ref. 4 and 5. 

The aim of this paper is to discuss another pathological case: the classi- 
cally integrable triaxial rotator model (see, for instance, Ref. 6). Incidentally, 
this model has been used very often in the description of the low-lying states 
of the even-even atomic nuclei.0 

The asymmetric top described by the rotor model is a classically inte- 
grable system, but an analytical formula for its quantum energy spectrum it 
not known. Nevertheless, numerical results can be obtained. By following 
the Landau approach, the Hamiltonian operator is split into 4 submatrices, 
corresponding to different symmetry classes. Each truncated submatrix is 
numerically diagonalized. Finally, the nearest-neighbor spacing distribution 
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P(s) and the spectral rigidity A3(L) are calculated. Surprisingly, the spectral 
statistics of energy levels do not follow the predictions of Poisson statistics. 
This result may suggest the presence of a hidden symmetry in the system but, 
to our knowledge, such a symmetry has not yet been identified. 

2 The Asymmetric Rotor Model 

Let us consider a system of coordinates with axes along the three principal 
axes of intertia of the top, and rotating with it. The classical Hamiltonian H 
of the top is given by 

H = \ {aJl + bJl + cJl) , (1) 

where J = {Ji, J2, J3) is the angular momentum of the rotation and a = 
b = I//2, c = I//3 are three parameters such that /i, I2 and I3 are 
the principal momenta of intertia of the top. The Hamiltonian is classically 
integrable and its action variables are precisely the three components Js, 
s — 1, 2, 3, of the angular momentum. 

The quantum Hamiltonian H is obtained by replacing the components 
of the angular momentum, in the classical expression of the energy, by the 
corresponding quantum operators Ji, J2 and J3. The commutation rules for 
the operators of the angular momentum components in the rotating system 
of coordinates are given by 

Jr'Js JsJr — ^rst Jt i (2) 

where e^st is the Ricci tensor and s, t = 1,2,3. Note that these commutation 
rules differ from those in the fixed system in the sign on the right-hand sidcH 
As usual, the two operators = Jl + J2+ Ji and J3 are simultaneously 
diagonalized on the basis of eigenstates | J, k) with integer eigenvalues J and 
k [k — —J,~J+ 1,...,J — !,•/), respectively. The non-zero matrix elements 
of the quantum Hamiltonian H in the basis | J, k) are given by 

(J, k\H\J, k) = ^{a + h){J{J + 1) - e) + ycfc2 , (3) 
(J, k\H\J, k + 2) ^{J,k + 2\H\J, k) = 

= ^[a-b)^{J -k)[J -k-l){J + k + !){.] + k + 2) . (4) 


The quantum Hamiltonian H has matrix elements only for transitions with 
/c — > fc or fc± 2. The absence of matrix elements for transitions between states 
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Table 1. Number of states in each submatrix of the asymmetrical top Hamiltonian for a 
fixed J. 
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J J 

2 2 


J odd 


J-1 
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.7+1 
2 


J+1 J+1 
2 2 



with even and odd k has the result that the matrix of degree 2 J + 1 is the 
direct product of two matrices of degrees J and J+1. One of these contains 
matrix elements for transitions between states with even fc, and the other 
contains those for transitions between states with odd fcfl 
It is useful to introduce a new basis, given by 

\J,k;S)^^{\J,k) + \J,-k)) , |J,0,5) = |J,0) , 

\J,k;A)^^{\J,k)-\J,^k)) , k^O. (5) 

By using this new basis, the total Hamiltonian matrix is decomposed in the 
direct product of 4 submatrices by considering the parity of the quantum 
number k: even (E) or odd (O), and the symmetry of the state: symmetric 
(S) or anti-symmetric (A). So the submatrices are labelled as follow: (E,S), 
(E,A), (0,S), (0,A). These are the classes of symmetry of the system. In 
Table 1 we show the dimension of each submatrix for a fixed J. 

The matrix elements of the Hamiltonian H in the new basis, with respect 
to the old basis, are given by 

{J,k,S\H\J,k,S) {J,k,A\H\J,k,A) ^ {J,k\H\J,k) , k^l (6) 

(J, 1, S\H\J, 1, S) = (J, 1\H\J, 1) + (J, l\H\J, -1) (7) 

(J, 1, A\H\J, 1, A) = (J, 1\H\J, 1) - (J, l\H\J, -1) (8) 

{J,k,S\H\J,k + 2,S) {J,k,A\H\J,k + 2,A) = {J,k\H\J,k + 2) , k^O 

(9) 

{J,0,S\H\J,2,S) ^V2{J,0\H\J,2) , k^O (10) 

We calculate the eigenvalues of each submatrix for different values of J using 
a fast implenifintation, in double precision, of the Lanczos algorithm with a 
LAPAC code.l 
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In Figure 1 we plot the density of levels p{E) of each submatrix of H and 
J — 1000. The results show that the density of levels is practically the same 
for the four classes. p{E) displays a high peak at the left-center of the energy 
interval and a long tail for large energy values. 



3 Spectral Statistics 



As previously discussed, according to the Berry- Tabor theoryliJ'LlJ, given a 
classical integrable Hamiltonian that, written in action variables Jr, satisfies 
the condition gj ^ 0, then, in the semiclassical limit, its spectral statis- 
tics should follow the Poisson statistics. Note that a system of linear harmonic 
oscillators, whose Hamiltonian is given hy H = lu ■ I, does not satisfy the pre- 
vious condition. In fact, a system of linear harmonic oscillators is integrable 
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Figure 2. Nearest neighbor spacing distribution P{s) of the four classes of symmetry with 
J = 1000. The dashed Hne is the Poisson prediction P{s) = exp (—s). Parameters: a = 1, 
b = V2, c = V5 and h = 1. 



but it docs not follow Poissonian statistics 

The triaxial rigid rotator is intcgrable and satisfies the previous Berry- 
Tabor condition. Thus, one expects that the spectral statistics of the quan- 
tized rigid rotator should be Poissonian. We shall show that is not the case. 

In general, various statistics may be used to show the local correlations 
of the energy levels but the most used spectral statistics are P{s) and A^lL). 
P{s) is the distribution of nearest-neighbor spacings Si = (i?t+i — Ei) of the 
unfolded levels Ei. It is obtained by accumulating the number of spacings 
that lie within the bin (s,^j±_As) and then normalizing P{s) to unit. As 
shown by Berry and Tabort2lO, for quantum systems whose classical analogs 
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Figure 3. Spectral rigidity A3(L) of the four classes of symmetry with J = 1000. The 
dashed line is the Poisson prediction A3(L) = L/15. Parameters: a = 1, b = \/2, c = Vs 
and h = 1. 



are integrable, P{s) is expected to follow the Poisson distribution 

P(s) = exp (-s) . (11) 

The statistic A3(L) is defined, for a fixed interval {—L/2, L/2), as the 
least-square deviation of the staircase function N{E) from the best straight 
line fitting it: 

1 /-^/^ 
A3(i) = - min / [N{E) - AE - BfdE , 

L A,B J_L/2 

where N{E) is the number of levels between E and zero for positive energy, 
between —E and zero for negative energy. The A3(L) statistic provides a 
measure of the degree of rigidity of the spectrum: for a given interval L, 
the smaller A3(L) is, the stronger is the rigidity, signiiying the long-range 
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correlations between levels. For this statistic the Poissonian prediction is 

AsiD ^ ^ . (12) 

It is useful to remember that Berry, on the basis of the Gutwiller semiclassical 
formula for the density of states, has shown that A3(i) deviates from the 
universal Poissonian predictions for large L: A^^p) should saturate to an 
asymptotic value performing damped oscillationsJla 

In Figure 2 the spectral statistic P{s) is plotted for the four submatrices of 
H and J — 1000. Note that the level spectrum is mapped into unfolded levels 
with quasi-uniform level density by using a standard procedure described in 
Ref. 13. As expected from the previous analysis of density of levels, P{s) is 
practically the same for the four classes of symmetry. Moreover, P{s) has a 
pathological behavior: a peak near s = 1 and nothing elsewhere. 

Compared to P{s), the spectral rigidity A3(L) is less pathological. As 
shown in Figure 3, A^iL) follows quite well the Poisson prediction A^{L) = 
L/15 for small L but for larger values of L it gets a constant mean value with 
fluctuations around this mean value. These fluctuaticms becomes very large 
by increasing L, in contrast with the Berry predictionllj. Note that we have 
repeated the caculations with other values of a, b and c but the results do not 
change. 

The behavior of the density of levels p{E) and of the spectral statistics 
P(s) and As^L) does not change by changing the matrix dimension, namely 
the quantum number J. In Figure 4 we plot the density of levels and the 
spectral statistics for J = 2000 and J = 4000. 

Conclusions 

The main conclusion of this paper is that the asymmetric rotor is, like the 
harmonic oscillator, another pathological case with respect to the classical- 
quantum correspondence between integrability and Poisson statistics. In our 
opinion, the pathology of the asymmetric rotor model is more interesting be- 
cause, unlike the harmonic oscillator, the asymmetric rotor satisfies the con- 
ditions of the Berry- Tabor theory. The presence of hidden symmetries could 
explain the pathological behavior of spectral statistics but such symmetries 
have not yet been identified. 

V.R.M. is greately indebted to Y. Alhassid, M.V. Berry, S. Graffi and 
V. Zelevinsky, for enlightening discussions. L.S. acknowledges M. Robnik for 
fruitful conversations. 
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Figure 4. Density of levels p{E), Nearest neighbor spacing distribution P{s) and spectral 
rigidity AsiL) of the (B, S) class of symmetry, with J = 2000 (top) and J = 4000 (bottom). 
Dashed lines are Poisson predictions. Parameters: a = 1, 6 = V2, c = V5 and h = 1. 
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